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We explore the constraints following from requiring the Area Law in the entanglement entropy in the 
context of loop quantum gravity. We find a unique solution to the single link wave-function in the 
large j limit, believed to be appropriate in the semi-classical limit. We then generalize our consid¬ 
erations to multi-link coherent states, and find that the area law is preserved very generically using 
our single link wave-function as a building block. Finally, we develop the framework that generates 
families of multi-link states that preserve the area law while avoiding macroscopic entanglement, 
the space-time analogue of “Schrodinger cat”. We note that these states, defined on a given set of 
graphs, are the ground states of some local Hamiltonian that can be constructed explicitly. This can 
potentially shed light on the construction of the appropriate Hamiltonian constraints in the LQG 
framework. 
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In recent years, the study of entanglement has shed new 
lights in the fields of quantum field theory, quantum sta¬ 
tistical mechanics and condensed matter theory. When¬ 
ever quantum fluctuations are important, entanglement 
plays an important role [1]. In the context of condensed 
matter physics, it has played a key role in classifying 
phases of matter, notably those with mass gaps, leading 
to the notion of topological orders [2-5] characterized 
by long-range entanglement [6-9]. A crucial quantifier 
of quantum entanglement is the entanglement entropy. 
For a wave-function that is the ground state of a local 
Hamiltonian, it is believed that the entanglement entropy 
should obey the so called area law [10]: in the limit of 
a large subregion A , its entanglement entropy is given 
by Sa = — T^Pa In Pa, where pa is the reduced density 
matrix of subsystem A , and Sa scales to leading order as 
the area of the boundary of A. 

From lessons that have been drawn from black hole 
entropy, which is proposed to be understood as an en¬ 
tanglement entropy between degrees of freedom within 
and without the horizon [11-13], and more recently evi¬ 
dence from the computation of entanglement entropy via 
the Ryu-Takayanagi formula [14] in the AdS/CFT cor¬ 
respondence, it appears that one very important feature 
of quantum gravity at least in the semi-classical limit, 
in which some classical background geometry can be de¬ 
fined, is that the leading contribution to the entangle¬ 
ment entropy associated to a sub-region A should also 
satisfy the area law [15]. 

Loop Quantum Gravity (LQG) [16-18] provides a 
perfect arena to test the requirement of area law, in 
which explicit calculations can be done. A long stand¬ 
ing and central question in the LQG is to understand 
the senri-classical limit. Thus far, senri-classical states 
have been proposed [19-22] in LQG, and studied in many 
contexts[23-40]. However, it is still by far an open ques¬ 
tion whether a semi-classical geometry has been undis¬ 


putedly recovered. The current paper is partly inspired 
by recent work in LQG that recovers the Bekenstein’s 
black hole entropy [41-50] and several attempts in com¬ 
puting entanglement entropy in terms of spin network 
states [51-55], which is intimately related to the area law 
of the entanglement entropy [15, 55]. 

In this letter, we inspect the semi-classical limit 
through the lens of entanglement entropy of a bounded 
sub-region in LQG. We find that imposing the area law 
across a single triangular surface, the simplest possible 
scenario, already leads to a very tight constraint on the 
possible wave-function, which admits an almost unique 
solution, one that acquires an interesting correction from 
what is envisaged in [15]. 

Ultimately the main goal of this paper is to show that 
entanglement imposes very stringent requirement to a 
theory of quantum gravity. We show how to use this 
guide in a constructive way. In particular, we impose 
that, together with the area-law, entanglement has to 
obey (i) SU(2) gauge invariance, and (ii) it must be mi¬ 
croscopic, implying that fluctuations in quantum space- 
times should be confined to microscopic scales. These 
are very natural requirements if Entanglement is phys¬ 
ically measurable. Gauge invariance means that entan¬ 
glement does result in stronger than classical correlations 
in actual physical observables [56], and microscopic en¬ 
tanglement means that we do not allow arbitrarily large 
Schrodinger cat-spacetimes [57]. These considerations, 
together with the single-link result, allows us to construct 
the many-body or multi-facet state of a general quantum 
geometry. As a bonus, we obtain a constructive method 
to build local Hamiltonians for given set of graphs, poten¬ 
tially important towards understanding the Hamiltonian 
constraints of LQG. 

Let us begin by setting the stage and notations. Our 
starting point is the SU(2) spin-network states in the 
graph-fixed kinematical Hilbert space Hr of LQG. An 
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SU(2) spin-network state is a triplet |«S) = \T,ji,i n ). T 
is the given proper graph with L oriented links and N 
nodes, ji € Z + /2 is an assignment of an SU(2) unitary 
irreducible representation to each link l, and i n is an as¬ 
signment of an SU(2) intertwiner to each node n. The 
spin-/ representation is isomorphic to a SL(2,C) repre¬ 
sentation ( 7 ji,ji) selected by the linear version of the 
simplicity constraints [58], namely K = 7 L, in which K 
and L are boost and rotation generators, respectively. 
The simplicity constraints are used in the spin foam for¬ 
malism [59-62] for recovering gravity amplitudes from 
BF(topological) amplitudes[63, 64] and for regaining the 
real connection [62, 65]. Spin-network states correspond 
to the so called twisted geometries [ 66 , 67], which de¬ 
scribe the geometries of three dimensional fuzzy discrete 
manifolds. Each IV-valent node corresponds to a iV-facet 
polyhedron, while each link is dual to the face of the 
polyhedron [ 68 ]. The areas of the faces are realized as 
the expectation value of the SU(2) Casimir operator act¬ 
ing on the spin-network states, i.e. the area of the face 
dual to link l is 


(Ai) = 8777 i 2 v yjji{ji + 1 ), ( 1 ) 

where 7 £ R + in the rest of the letter, and i p is the 
Planck length. 



Figure 1. A face Ai dual to link l. 


Given a state |'F) in Hr and cut the graph T into two re¬ 
gions A and A, the entanglement entropy Sa between 
these two regions is defined as the von Neumann en¬ 
tropy of the reduced density matrix pa = Tr^ |'I') {T |. 
Tr A is the partial trace over the states (i.e. the inter¬ 
twiners’ and spins’ degrees of freedom) inside A. To re¬ 
cover locally a Minkowskian vacuum that respects the 
Bisognano-Wiclnnann theorem, it is proposed that the 
reduced density matrix pA should be proportional to 
the exponentiation of the boost operator K , leading to 
the following proposal for a single link state [15, 69] 

14/) = Ylj a j exp(— jK) mn \j, to, n). A careful check how¬ 
ever suggests that it does not satisfy the area law [70] 
in the large j limit, which is another crucial ingredient 
of the semi-classical limit, opposing naive expectation. 
Indeed, one can readily check that the normalization of 
the reduced density matrix always cancel the area term 
in this limit. Surprisingly, in the large j regime, there is 
a natural and arguably unique solution that recovers the 
area law across each link. For a single link we can conve¬ 
niently choose a state in the diagonalized basis with fixed 


spin j such that the reduced density matrix satisfies 
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E P( m )\A rn )(A m \’ P( m ) 

m——j 


f M 


The function f(m) is assumed not to carry explicit j 
dependence. We will also replace the sum by an integral 
J2m=-j f-j dji assuming that the distribution p(m) 
is negligible at small values of |?ti| in the large spin limit, 
an expectation that is naturally born out in the original 
proposal of Ref. [69]. Requiring that the entropy admits 
an area law and taking into account ( 1 ), we impose Sa = 
cj + • • •, the derivative of which w.r.t. j is thus a constant 
c > 0. This leads to a constraint on f(m) [71]: 

+ cj -1 + In p(j) = 0 . ( 2 ) 

PW 

Eq. (2) cannot be solved exactly in complete generality. 
However, anticipating that we are considering the large 
spin limit, and motivated by the original proposal, we 
consider the limit in which c/p(j) -C \cj — 1|. In this 
limit, we have p(j) = exp(l — cj), which finally gives 


/(to) = exp 


—CTO 


exp(l — cto) 
c 


(3) 


One can readily check that c/p(j ) <C | cj — 1| is satisfied 
in the large j limit if we require that c scales as 

0 < j- 1 « c« j~i < 1. (4) 


In this limit the leading term in the entanglement en¬ 
tropy is indeed linear in cj, recovering the long sought 
area law. The function p(m) also naturally suppresses 
contributions at small values of to in this limit, justifying 
the approximation that replaces the sum by the integral. 
It is very interesting that the state that we have found is 
very similar to the original proposal, up to an exponen¬ 
tially suppressed factor, which is negligible where /(to) 
actually contributes, but suppresses regions which would 
otherwise have contributed in the original proposal. The 
constant c also emerges naturally and plays precisely the 
same role as the Barbero-Immirzi parameter 7 , and (4) 
is exactly the semi-classical large j regime in the covari¬ 
ant formalism of LQG [31]. Recall that the boost oper¬ 
ator K z was related to L z via the simplicity constraint 
K = 7 L, which followed from an action quadratic in 
the tetrad. One might be tempted to interpret the ex¬ 
tra exponential correction we find here as a non-linear 
correction arising from quantum effects. Combining our 
proposal with the well known construction of multiple 
link coherent spin-network state [19-22], our new state 
14') in Hr is found to be 

= ^ E (5) 

lizl.Tr,} l,n 
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In (5) we defined: A-’A j^ = dj t e~^ l ~ Jl ^ tl , with dj t = 
2ji + 1 ; Ji, a given spin for the link l , and t;, similar 
to a heat kernel time; the Wigner matrix of the group 
elements that encodes the correction found in (3), i.e. 


p.h = pit 



-K~/Lz+icf>iL z - 


exp(l — 2 n'yL z ) j.\ 

9i) , 


( 6 ) 


in which • •) is defined as (j, m\ • • ■ \j, n ); <fn is a 

phase; gi and gi are SU(2) elements; L z is the su(2) gen¬ 
erator in the z axis; ly is the normalized intertwiner with 
Vi ■ Vi* = 6u>. Note that this is a state that is factoriz- 
able into wave-functions of individual links up to gauge 
constraints. 

We can extend our construction to the more compli¬ 
cated case of one node, with many links protruding from 
it. For simplicity, the graph T used in the following cal¬ 
culation is the graph with L out-pointing links attaching 
to only one node. 



Figure 2. A graph with L links and one node. 


We cut the graph into two regions as in Fig.2. Region 
A contains the node. The reduced density matrix pa is 
recovered by tracing out the link source states (or the 
intertwiner degrees of freedom) in |’F), which reads 


i*) = S E <"> IFi 
N m Jsum Vi 




x 

{ni} l 


Furthermore, pa can be rewritten in the block-diagonal 
form p A = E {jl }P{ji} = where each 

block is characterized by the configuration of the {;/)} 
eigenvalues of the boundary links. Crucially, such a block 
factorization follows from the SU(2) invariance imposed 
at the node and the properties of the intertwiners. We 
can write P{j t \, the norm of the {jz}’s block P{j t }, as 


p m = jp J sv ^ ( 8 ) 

In ( 8 ) we use that: f(L z ) is given by (3) with c = 27ry; 
Xj, (M) is the trace over matrix M; p^^ is the reduced 
density matrix with fixed {j{\ 


P{ji} [ dhJ[ pJl ■ D{h) • F* t, (9) 

Z {ji} J SU(2) , 


where F Jl is given by ( 6 ); Zypi is the normalization such 
that Trp = 1. The integration over SU(2) element h fol¬ 
lows from the summation of intertwiner, which is propor¬ 
tional to an integration over Wigner matrices on SU(2) 

E, 4 P ‘ } ^ Pl} « fsu( 2 ) dh Ut D p iv [{ h )■ The normaliza¬ 
tion has been absorbed into A f. Then the entanglement 
entropy Sa = —Tr pA In pA evaluates to 

S A = - J2 P {v} Tt P{ji} ln P{n} - P {v} lnP {M ■ 
{ni m 

Note that: i) P{j t \ can be understood as the probabil¬ 
ity density of detecting the given boundary configuration 
{ji}, and that ^{j,} P {n} = 1 > b) it is a subtle issue to 
define the entanglement entropy in a gauge theory [72, 73] 
(moving from the proposal in [15], the resultant Sa is in 
fact a SU(2) gauge invariant quantity which only has de¬ 
pendence on the eigenvalues of the Casimir of the bound¬ 
ary links, a choice equivalent to the proposal in [74]); iii) 
the entanglement entropy with a unique boundary link- 
configuration is given in the first term of ( 10 ), and it will 
be proved later to be the area term in the large spin limit. 

To calculate Sa, we first give an explicit expression 
for PtjA and Tr/b^i ln/Sr^i. As such, we introduce the 
coherent states on SU(2) [21], namely D J m -(n) \j,m) = 

n\j,j) = n |t)' gl2j , with n £ SU(2). Imposing the SU(2) 
coherent states resolution of identity into each \ji i n (8), 
P {ji} becomes, 


p m = jp I dhY[dm a\\ u 


(— e ) ? 


Jl) Jl p Pi'-( 27T7)p< 


where the “action” S p is given by 


L 

S P = Y, Vi In (tl njme-^+^gjh m |t) • (10) 

i 


In order to study the large spin behavior of P{j t }, we 
rescale ji —► Xki and assume A 1. It is convenient to 
perform the asymptotic expansion to get the major con¬ 
tribution of P{j t }- The solutions of the equations of mo¬ 
tion 6 ni S p = ShSp = 0 control the senri-classical behavior 
of PtjA- For S ni S p = 0, the only possible solutions [75] 
are ni = gi, h = 1. Note that ShS p = 0 is equivalent to 
the Gauss constraint. Then P{j t }, to leading order of the 
asymptotic expansion and in the limit of (4), becomes 


T \n} = Jp\{ A 0 „ J ,)^ ex P(- 1 - 


2 ji 


a l-27T7 j 


2ttj 


( 11 ) 


In order to estimate the term Tr/3| 71 | In P { :n \, it is 
easier to use the replica trick to compute S n , the 
Renyi entropy of order n [76]. From Trp7., = 

nfiz n /(^)e ( ”- 1)(1 - 2 ^ !) and taking the n — > 1 limit, 
we get -Tr p {jl} In p {jl} = 2ivyj t - In (jlf jl) ■ Fi¬ 

nally, we have 
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Sa =X^(2tt7 ji) - yin J - ( lnP {oi }) > ( 12 ) 

where we have denoted (■■■} = (■ ■ • )P{j}- The first 

term is the area A of region A’s boundary. The loga¬ 
rithmic correction is coming from both the second and 
third terms. When ti = J^ k , the logarithmic correction 
is (—k + 2) In J;. So in the semi-classical limit of LQG 
(4), the entanglement entropy is 

^ = (13) 

where the “chemical potential” term now carries only 1 / j 
suppressed contributions. We would like to generalize 
our construction to multi-node states. Let us start with a 
simple state that still satisfies the area law by considering 
just the tensor product of single node states over all the 
nodes of the graph. 

\^n) = Y,{j,},{i n } Hoi }{*«} Yin l^n), (14) 

I'l'n) = \{jlen},in) , c {ji}{i n } = Jf II/, n ^(t;, Ji)^ ' Vi n ■ 

The links that are completely enclosed within the re¬ 
gion A would not contribute to the entanglement entropy, 
which will instead scale with the number of links crossing 
the boundary dA yielding the area law, A = YYiedAi^l) ■ 
From the many-body point of view, this state is triv¬ 
ial, as there are no quantum fluctuations other than the 
single-body ones. A generic superposition of such kind 
of states, though, will first of all not respect the area 
law, and will contain arbitrarily macroscopic superpo¬ 
sitions. We need to deform the trivial state in a way 
that some entanglement and quantum fluctuations are 
produced, but without violating the constraints we set 
up. This can be obtained by the technique of quasi- 
adiabatic continuation introduced in [77]. This technique 
allows to continuously deform a quantum state that is 
the ground state of a local Hamiltonian in order that it 
is still the ground state of a local Hamiltonian. Let H( A) 
be a smooth family of local Hamiltonians parametrized 
by A € [0,1] such that |fEw) = l'VE'(O)) is the ground state 
of H(A = 0). Here, local means that H( A) is the sum of 
local operators, H{A) = YZ X Txcr(A), where X denotes 
a subset of the graph T and Tx{ A) means that this op¬ 
erator has only support on X. Also assume that H{ A) 
has a finite gap A E between ground and first excited 
state for all A’s. One way to think of H{ A) is as the 
perturbation of the initial Hamiltonian H( 0), namely, 
H( A) = 17(0) + X^2 x Kx and K x are local operators 
with support on X. Then, following [78], one can de¬ 
fine a unitary operator U{ A) = TexpJ— i f* H(s)dsj 
with T the time ordering operator, and 7f(s) = 

i J dtF(t) exp(iH(s)t)(d s H(s)) exp(— i.H(s)t), with F(t) 
an appropriate fast decaying smooth function (see [78] 


for details). The unitary operator U{ A) so defined in¬ 
terpolates among ground states of 17(A) and it is called 
adiabatic continuation [77]. The ground state of 17(A), 
l’F(A)), can be written as the adiabatic continuation of 
[4/(0)), namely l'P(A)) = U{ A) [^(O)). Now, the adia¬ 
batic continuation preserves the area law [10]. In this 
way, we can always deform |'L( 0 )) in a way to obtain a 
new state | (A)) such that its parent Hamiltonian 17 (A) 
is non-integrable. 

We now show that by adiabatic continuation we also 
preserve the physicality of entanglement as being micro¬ 
scopic. Adiabatic continuation preservers macroscopic 
entanglement. So, if one continues a state that is not 
a Schrodingers’ cat, one will not obtain a Schrodingers’ 
cat. On the other hand, continuing a Schrodingers’ cat 
will still yield a macroscopically entangled state. To show 
this, we apply a result obtained in [57]. The superse¬ 
lection rule is imposed by stating that states with non 
vanishing mutual information between two distant 
macroscopic regions must be ruled out. The mutual in¬ 
formation is defined asF(A\B) := S(A) + S(B) — S(AB). 
Here, S is chosen to be the 2—Renyi entropy S -2 = 
— logTrp^, namely, the logarithm of the purity instead 
the von Neumann entropy. The technical reason for mak¬ 
ing this choice is that one can prove that I oo is preserved 
by quasi adiabatic continuation. The physical reason, 
is that this quantity is a physical observable (unlike the 
Von Neumann Entropy) [79-81], but still captures all the 
entanglement properties that are important in quantum 
many-body theory [82, 83]. As a bonus, starting with 
a trivial non-interacting Hamiltonian that is the sum of 
commuting local terms in the tensor product structure of 
Eq.(14), one can adiabatically continue the very Hamil¬ 
tonian by means of any other Hamiltonian with a gap, 
which is the sum of local operators with couplings de¬ 
pending on some parameter A. This is thus a systematic 
construction of local Hamiltonians whose ground states 
satisfy the desired properties we have set forth. 

To conclude, we have shown that entanglement can 
provide stringent guiding principles in selecting states 
with sensible semiclassical limits in LQG. We obtain 
these results by imposing that entanglement must be 
physical (gauge invariant and not macroscopic) by using 
methods from quantum information and quantum many- 
body theory. In perspective, the framework here devel¬ 
oped will allow for the study of notions like quantum 
order or thermalization in a closed quantum system in 
the context of quantum gravity. 
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